Optical selection rules of graphene nanoribbons 
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Optical selection rules for one-dimensional graphene nanoribbons are analytically studied and 
clarified based on the tight-binding model. A theoretical explanation, through analyzing the velocity 
matrix elements and the features of wavefunctions, can account for the selection rules, which depend 
on the edge structure of nanoribbon, namely armchair or zigzag edges. The selection rule of armchair 
nanoribbons is A J = J" — J" — 0, and the optical transitions occur from the conduction to valence 
subbands of the same index. Such a selection rule originates in the relationships between two 
sublattices and between conduction and valence subbands. On the other hand, zigzag nanoribbons 
exhibit the selection rule |AJ| — odd, which results from the alternatively changing symmetry 
property as the subband index increases. An efficiently theoretical prediction on transition energies 
is obtained with the application of selection rules. Furthermore, the energies of band edge states 
become experimentally attainable via optical measurements. 

PACS numbers: 78.67.-n, 73.22-f, 81.05.ue 
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INTRODUCTION 

The discovery of graphene [l|, Q has motivated many 
studies because of its interesting physical properties 
and possible applications. Graphene, a two-dimensional 
carbon-based material 0-0]' is a gapless semiconductor 
with linear energy dispersion near the Dirac point [^, 0| ■ 
This cone-like energy spectrum results in particular elec- 
tronic and transport properties, such as massless Dirac 
fermions which travel at a extremely high constant 
speed [1, finite conductivity at zero charge-carrier 
concentration [l^, 11|, and unusual quantum Hall ef- 
fect [H, [l^l- However, the gapless feature of graphene 
limits its applications in electronic devices. In order 
to induce a tunable band gap for a wide variety of ap- 
plications, one-dimensional (ID) graphene nanoribbons 
(GNRs) are proposed, which can be fabricated using 
a nanowire etch mask 14, U^, chemical vapor deposi- 



tion 

(Sin, 

lithographic patter ning of graphene 
and unzipping of nanotube 2ll-l24|. The electronic and 



transport properties are highly related to edge structure 
and ribbon width. There are two basic types of GNRs, 
armchair (AGNR) and zigzag (ZGNR) ones, according 
to their edge structures along the longitudinal direction. 
ZGNR is gapless because of the peculiar partial flat bands 
at the Fermi level {Ep — 0), while AGNR can be ei- 
ther semiconducting or gapless depending on the ribbon 
width. Owing to these special electronic properties, GNR 
has triggered many studies, e.g., geometrical and energy 
band structur es 12 5^, 26] , as well as magnetic IIQ, i27.-30] , 
transport 3ll433| , and optical properties 3J, |35| . 

Optical measurements are very useful in understanding 
the electronic structure of a material. Unlike the feature- 
less absorption spectrum of graphene at low energy [s^ , 
the optical spectra of GNRs exhibit many sharp peaks 
relating to van Hove singularities 37|, |38| in the density 



of electronic states. These peaks satisfy certain selec- 
tion rules, which depend on the edge structure. Absorp- 
tion selection rules of GNRs, which crucially decide the 
possible transitions, have been studied by Lin et. al. 
firstly 35]. Although these rules are described, the ori- 
gin is still not verified. In this work, the selection rules 
are analytically studied and a theoretical explanation is 
given. 



TIGHT-BINDING MODEL AND ELECTRONIC 
PROPERTIES 

The geometric structures of AGNR and ZGNR are 
shown in Fig. [1] Atoms A (dots) and B (circles) are 
the nearest neighbors to each other with the C-C bond 
length b = 1.42 A. The ribbon width is characterized 
by the number (Ny) of dimer lines (armchair or zigzag 
lines) along the longitudinal direction (x). The peri- 
odic length is h = 36 (^ = ^/36 = a) for AGNRs 
(ZGNRs). The first Brillouin zone is confined in the re- 
gion —n/Ix < kx < 1^ I Ix with 2Ny carbon atoms in a 
primitive unit cell. The Hamiltonian equation of the sys- 
tem is H\^) = E\^). With the framework of the tight- 
binding model considering only the nearest-neighbor in- 
teractions, the Hamiltonian operator is H = j, 7oc|c/', 

wliere 70 is the hopping integral, cj (c/') is the creation 
(annihilation) operator on the site I {I'), and the on-site 
energies are set to zero. The Bloch wavefunction, i.e., the 
superposition of the 2Ny tight-binding functions \(l)i)'s 
with the corresponding site amplitudes C/'s, is expressed 
as 
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FIG. 1. Energy dispersions and geometric structures for 
Ny = 69 AGNR and ZGNR. The conduction (valence) sub- 
bands of the corresponding indices J'^ (J") are depicted in red 
(blue) color. The dashed-line rectangles represent primitive 
unit cells. The numbers of mth dimer line are on the right 
side of the hexagonal lattice. The notations in the rectangles 
denote the A and B atoms on the mth dimer line. 



For ZGNRs, the Hamiltonian elements are 



270 cos(fc2,a/2), + / is odd, 

Hiji = ■{70, ii I' = I + 1, I is even, (3) 

0, otherwise. 



where 70 = 2.598 eV [39]. The eigenvalues £''"s and 
eigenstates |^''')'s are obtained by diagonalizing the 
Hamiltonian matrix. The superscript h can be either 
c or V representing the conduction or valence subband. 

The energy dispersions of GNRs are symmetric with 
respect to fc^; = 0, and therefore only the positive ones 
are discussed. For the AGNR of Ny — 69, many ID 
parabolic subbands due to the finite width of GNR are 
symmetric about the Fermi level [Fig.[lja)]. The bottoms 
and tops of these subbands are at fc^; =0, where many 
band-edge states exist and optical transitions may occur. 
The energy gap is equal to 0.05 70 for this semiconduct- 
ing nanoribbon. The subband index J'^'" {— 1, 2, 3, ...) is 
decided by the minimum energy spacing between the sub- 
bands and the Fermi level. For the ZGNR of Ny ^ 69, the 
parabolic subbands are also symmetric about Ep = 0, 
and their band edges are located at = 27r/3 [Fig.[IJb)]. 
Furthermore, there are degenerate partial flat bands ly- 
ing on the Fermi level due to the zigzag edge structure, 
and the dispersionless region extends from 27r/3 to tt. 
Hence, it is a gapless semiconductor. 



where A^n (Bm) is the site amplitude and \am) (|^m)) is 
the tight-binding function associated with the periodic A 
(B) atoms. 

The Hamiltonian H in the subspace spanned by the 
tight-binding functions is a 2A^j^ x 2Ny Hermitian matrix. 
For AGNRs, the upper triangular elements of H are 



1,1' 



70 exp[ifc:r (-&)], 
7oexp[ifc:r(-V2)], 
7oexp[ifc^(&/2)], 
0, 



HI' = 1 + 1, I is odd, 

ii I' = I + 1, I is even, 

if = / + 3, lis odd, 
otherwise. 



(2) 



FEATURES OF WAVEFUNCTIONS IN REAL 
SPACE 



Wavefunctions give information on the charge density 
distribution, and their features are very important in un- 
derstanding optical selection rules. They are strongly 
dependent on the edge structure of nanoribbon. For AG- 
NRs, the wavefunction is a combination of six subenve- 
lopc functions 
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(4) 



where A{B)i^„^'s, A{B)^j^j^i^s and A{B)i^,^^2'^: ^-re the (3to + l)th and {3m + 2)th armchair lines with m being 
site amplitudes of sublattice A (B) located at the (3TO)th, 
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FIG. 2. Wavefunctions of subband index J'^ — 1-3 (red dots) 
and J" = 1 and 2 (blue dots) at fc^ = for Ny = 69 AGNR. 



FIG. 3. Wavefunctions of subband index J'^ = 1-3 (red dots) 
and ,r = 1 and 2 (blue dots) at = 27r/3 for A^^ = 69 
ZGNR. 



an integer. The wavefunctions of Ny = 69 AGNR at 
kx = 0, where the optical transition occurs, are shown in 
Fig. [2] They have an oscillating pattern, and the num- 
ber of nodes (zero points) is given by — 1. The square 
of wavefunction is the charge density distribution. The 
local maxima and minima of wavefunctions are of high 
charge density and the nodes are zero. As the state en- 
ergy increases, the oscillations become more severe. 
More importantly, each of them exhibits two aspects of 
features. For a certain wavefunction, there is a relation 
between the subenvelope functions of sublattices A and 
B, that is, the site amplitudes A'^{J^) and B^{J^) abide 
by the relations 

At{f^)=±Bt{J% (5) 

For example, the site amplitudes of sublattice A of = 1 
subband, A^(J'= = 1) [Figs. [i;m)-12i;o)] , is the negative 
oiB^{.r = 1) [Figs.[2i;p)-[2i;r)]. For the conduction and 
valence wavefunctions of the same index, J'^ = = J, 
their subenvelope functions have the relations 

A^JJ) = ±Al^{J). (6) 

For instance, A'^{J — 2) — —A^{J — 2) [compare 
Figs, [2^g)^2^i) with Figs.[2i;y)-I2i;aa)]. The relations be- 
tween the subenvelope functions of the same subband 
index are either in phase or out of phase. 



For ZGNRs, the wavefunction is decomposed as 

odd even odd even 

(7) 

where A^"s, B^'s, A^'s and B^'s are the site amplitudes 
of sublattices A and B on odd and even zigzag lines. 
[ao)'s, [60) 's, [ae)'s, and |6e)'s are the corresponding tight- 
binding functions. The wavefunctions of Ny — 69 ZGNR 
at kx = 27r/3, as shown in Fig.jSj reveal profiles different 
from armchair ones. The subenvelope functions of the 
conduction and valence subbands follow the relations: 

where the indices {Ny + 1 — m) and m are correlated 
and the sign factor (—1)"' [(— l)'' ] describes whether 
the wavefunction is symmetric (positive sign) or anti- 
symmetric (negative sign). For example, the subband 
of J'^ = 1 is symmetric [compare Figs. |31i), with 
Figs. [3lk), [Sjl)] and anti-symmetric for J" = 1 [compare 
Figs. [SUm), [3fn) with Figs. [3lo), [Hp)]. It is notewor- 
thy that the symmetry property is associated with the 
subband index. Considering the conduction subbands, 
the J'^ = 1, 2, and 3 wavefunctions are symmetric, anti- 
symmetric and symmetric, respectively [Figs. |31a)-(3Ul)]. 
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In other words, the symmetry property changes alterna- 
tively from symmetry to anti-symmetry with the incre- 
ment of subband index. Moreover, the results of this 
study show that the characteristics of the subenvelope 
functions A'^{J'^) and B^{J'^) are very sensitive to Ny 
being odd or even. If Ny is even, the relations become 



which differ from Eq. ([8]) by a minus sign. However, 
the alternatively changing symmetry property remains 



whether Ny is odd or even. 



ABSORPTION SPECTRA 



AW) = i- 



1)"' B'f^y + l-niiJ")^ 



(9) 



When GNR is under the influence of an electromag- 
netic field with polarization E||x, the absorption spec- 
trum A{uj) for direct photon transitions (Afc^; = 0) from 
state \'^''{k^,j'')) to state |*''' (fc^;, j'*')) is given by 



A{uj) « E / 



h,h',J'^,J'^' ' 



1st B.Z. 



dkx 
"27 



(*'''(fc„j''')|^|*^fc,,j'')) 



X Ini ■ 



f[E\k,,J^)]~f[E'^'{k,,jf^')] 
Eh'{kx, Jh') - Ef^ik,, J>')-Lo- iT 



, (10) 



where vrie is the effective mass of electron, P is the mo- 
mentum operator, j'')] is the Fermi-Dirac dis- 
tribution function, and T is the broadening factor. The 
low-energy absorption spectra for Ny = 69 AGNR and 
ZGNR at zero temperature are illustrated in Fig. SI The 
spectra show a lot of ID peaks caused by the inter-band 
transitions of ID subbands, and they are greatly affected 
by the geometric structure. 

For absorption spectrum of the Ny — 69 AGNR 
[Fig. mja)], the first peak, located at wi ~ 0.055 70, 
is numerically identified as the excitation from — 1 
to = 1 subbands. The second (aj2 = 0.105 70) and 
the third (0-13 — 0.205 70) ones are the transitions from 
J" = 2 to = 2 subbands and from J'' = 3 to J^^ = 3 
subbands, respectively. The allowed transitions originate 
in the subbands of the same indices, i.e., the selection 
rule is A J = J'^ — = 0. The transition energy is twice 
the band-edge energy {ujj — 2E'^ff'^ — 2\E^ff^\, where 
J'^ = J"" — J) because the band structure is symmetric 
about Ep = 0. Moreover, the peak arises higher when 
the frequency increases as a consequence of the decrease 
of subbands' curvatures. 

For ZGNRs, the absorption peaks may be classified 
into principal peaks and subpeaks according to their 
heights, as shown in Fig. IH^b). The allowed transitions 
obey the selection rule |AJ| = odd, which is different 
from the armchair case, AJ = 0. The principal peaks 
become higher with the increasing frequency because of 
the decrease of bands' curvatures and the increase of ex- 
citation channels. Due to the symmetry of energy spec- 
trum, the absorption spectrum of negative AJ's is the 
same as that of positive ones. Therefore, only the pos- 
itive AJ cases (transitions from J" to = J" -|- AJ 
subband of energy ui — Eji^^^j — E'^fj'^) are discussed. 
The principal peaks, Wj , are contributed by transitions 
between subbands satisfying the selection rule AJ = 1 



{J" J" : I 2), AJ = I (J" ^ J^ : 2 3), 
AJ = 1,3 (J^ ^ J^ : 3 ^ 4, 2 ^ 5,), A J = 1,3,5 
(J^ ^ J'^ : 4 5, 3 ^ 6, 2 ^ 7), and so on. The 
corresponding numbers of excitation channels are 2, 2, 4, 
and 6, i.e., the double of the allowed positive A J number. 
On the other hand, there exists one subpeak between two 
adjacent principal peaks. These subpeaks are indicated 
by wij, where the subscript IJ denotes the transition 
from the J" = I valance subband to the J"^ conduction 
subband. The subpeak W14, for instance, results from 
the transition between the J"" = 1 and J"^ = 4 subbands 
and follows the optical selection rule A J = 3. Since each 
subpeak possesses two excitation channels, the increase 
of subpeak height is not as fast as the principal one. It 
is remarkable that only the first principal peak (wf ) and 
all the subpeaks (wij) are from the transitions related 
to the J='^ = I subbands {Ejc.^i = Ep = 0). The 
peak is located at E^ff^^^ ^-^d the subpeaks wij are posi- 
tioned at Ejt^J^j, which is the energy difference between 
the subband edge and the Fermi level. 

OPTICAL SELECTION RULES OF GNRS 

In addition to the above-discussed features of wave- 
functions, the optical selection rules are also deter- 
mined by the velocity matrix elements [M'''* (k^) = 
(^''''(fc^, J''')|(E • P/me)|*''(fc:r, J'*))]. To evaluate M, 
the gradient approximation (E • P/rrie — V^iJ) is em- 
ployed [i^l, so that 

M'''''(fc.) = J2 C^(fc.)Q'?'(fc.)^^||^. (11) 

i,i'=i ^ 

At zero temperature, only the inter-band transitions 
from the valence to conduction subbands are valid. To 
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FIG. 4. Absorption spectra for Afy = 69 AGNR and ZGNR. 



investigate the optical properties, we focus on = 0, 
where the excitation occurs for AGNRs. Substituting 
the sublattices B with A in Eq. ([TT]) by the relations 
B-^ = s< (s = ±1) and B^^ = tA^^ {t = ±1) from 
Eq. ©, M is 



m—1 



m—1 

E 

m=2 



2 

670 



- tAZ{nAi^+,{r) + sAZ{nAi^+,{r) 

- tAZ{r)Ai,_,{r) + sAZ{nAi,_,{r) 

(12) 



The transfer relations between Am±i and Am, obtained 
from = E\'^') at fc^^ = 0, are 

A2 = A±Ai, 

+ = A± A™, if m = 2, 3, 4, A^j, - 1, 

An^^i = A±Ajv„, 

(13) 

where A± = ±(i?/7o) — 1 for Am = ii?™- After using 



Eq. p3)) to replace Am±i with A„i, M becomes 



^70 E 



sA-(J^)A;^(J")+tA-(J=)A:;(J'') 



-^tA';:ir)Ai,{r) + ^sAZ{r)A"m{,r) 



m—1 

where A^ = s{E/-iq) - 1 for A%^ = sB,'^^. Through the 
relation between the conduction and valence subbands 
ylj^j — uAm {u = ±1) from Eq. ^ and the normal- 
ization of the wavefunction '^Y.mLi A^{J'')A^^{J'') ^ 
(^J^.J") the velocity matrix element at fc^^ = is 



M'"=(0) = u{t-s){l 



±(1 - ^)67o^J-,J-, 



(15) 



where A± — ±(i?/7o) — 1 and u{t — s) is cither or ±2. 
Eq. (US]) achieves the selection rule (A J = 0) for AGNRs 
and possible excitation channels in absorption spectra. 

For ZGNRs, the velocity matrix element in Eq. pT|) at 
kx = 27r/3 becomes 



A';:ir)Bi{r) + B'^ir)Ai,ir) 



(16) 



where h' — — 070 sin(7ra/3). For odd Ny, we use Eq. ^ 
to convert the second term in Eq. ([T5]) . which leads to 



m — 1 

+(-1) 



A'^^ij^Bur) 



.(17) 



Then renaming the index of A^j^ + 1 — m as m gives 

Ny 

[1 + i-iy'+'^+'w ^ A';::{jnBiir). m 

m—1 

As for even Ny, the second term in Eq. (flT)) can also be 
converted by Eq. ([9]), and it yields the same result as 
Eq. (dll). Hence, M^'^ for ZGNR at = 27r/3 is 

M-r — ) = / 2/^' (./^)^m(^''), if A./ - odd, 

^ 3 ^ \ 0, if AJ = even, 

(19) 

where J2rn=i ^miJ'')BmiJ^) a non-zero value, and 
the selection rule (A J = odd) for ZGNRs is obtained. 

By using the selection rules, the transition energies 
(peak positions) can be efficiently obtained from the en- 
ergy dispersion. For AGNRs, the dependence of the first 
five predicted transition energies on the ribbon width is 
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ujie = Eji^^Q, and uj^ — E'^ffji^ — Eji^^^. The first sub- 
peak and the second principal peak are close and will 
merge for sufhciently large width, so do the second sub- 
peak and the third principal peak. In short, the peak 
frequencies can be predicted by the combination of band 
structures and selection rules. This is an efficient way to 
obtain wide-range information on the transition energies 
without extrapolation. 

It is worth mentioning that the exact energies of band- 
edge states for ZGNRs can be specified by the exper- 
imentally measured transition energy uj'y'. According 
to the optical selection rule, the measured peak positions 



P,exp 



exp 



14 , and Wig 

edge rpedge 



can be applied to specify the band- 
edge states Ejt^J^i ^j^=4j Ej^^J^, respectively. They 
are the energy differences between the J'' — even sub- 
band edges and the Fermi level. Moreover, the band-edge 
state energies of the J"^ = odd subbands are obtainable, 

r 1 TT^edne P,exp i rr^edne i P,exp P,exp 

lor example, iijei^3 = ^2 |^ji.=2l~^2 '^i ■ 
Similar backtracking strategy is valid for AGNRs based 
on the band symmetry about the Fermi energy. The 
band-edge state energy, which is the half of the exper- 
imentally measured peak position, is given by E'jt^'^ = 



/2, where J'^ = = J. Hence, the selec- 



l^edgei _ ^ exp 

tion rule provides a way to overcome the disadvantage of 
optical measurement, that is, it only yields the informa- 
tion about the energy difference between two subbands. 



FIG. 5. Tlie first five consecutive transition energies witii 
respect to the ribbon width Ny for AGNR and ZGNR. The 
notations •, o and ■ correspond to the AGNRs of Ny — 3m, 
3m -|- 1, and 3m -I- 2, respectively. 



shown in Fig. El^a). For a certain width Ny, the Jth 
transition energy is equal to u!j = 2E'^ff^ , i.e., twice 
of the band-edge state energy of the J'^th subband. For 
narrower GNRs, there are three groups of transition ener- 

1 

-2 



gies. The wj's of Ny = 3to (full circles) and Ny = 3m- 
(open circles) groups are close. While for the Ny = 3m- 
(squares) group, the first and second wj's are close and 
sandwiched between the second and third wj's of the 
Ny ^ 3m -I- 2 nanoribbons. For wider GNRs, these 
transition energies make red-shift and merge together. 
The GNRs may therefore be sorted into two groups of 
Ny = 3m -I- 2 and Ny ^ 3m -|- 2. The former is gapless 
due to the linear bands crossing at the Fermi level, while 
the latter is semiconducting with a band gap energy cor- 
responding to the first transition energy. 

The first five ribbon-width-dependent transition ener- 
gies of ZGNRs are presented in Fig.ISjb). They are asso- 
ciated with the first principal peak wf , the first subpeak 
wi4, the second principal peak wl', the second subpeak 
wie, and the third principal peak u)^ , respectively [see 
Fig. m^b)]. The predicted corresponding transition ener- 



CONCLUSION 

We employ the tight-binding model to study the ab- 
sorption spectra of GNRs, and this is for the first time 
that the corresponding optical selection rules are ana- 
lytically specified. The main results of this work are 
stated as follows. First, the optical transition channels of 
absorption spectra for AGNRs and ZGNRs are exactly 
identified. Then, the characteristics of absorption peaks, 
such as positions and heights, related to the energy dis- 
persions are discussed in detail. Most importantly, this 
work provides a theoretical explanation on the optical 
selection rules through analyzing the velocity matrix el- 
ements and the wavefunction features. The optical se- 
lection rules depend on the edge structure, armchair or 
zigzag edges. The selection rule of AGNRs (AJ ~ 0) 
originates from two aspects: the relation between sub- 
lattices A and B within a certain state and the relation 
between the conduction and valence subenvelope func- 
tions of the same index. While for ZGNRs, the selec- 
tion rule (I A J| — odd) owes to the alternatively changing 
symmetry property with the increasing index. Accord- 
ing to the selection rules, we can efficiently predict the 
peak positions and identify the corresponding transition 
channels. Furthermore, this gives a way to overcome the 
limitation of optical measurements. Hence, the energies 
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of band-edge states can be exactly obtained from optical 
experiments. 
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